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Abstract 

On the basis of the Edward-Kornfeld formulation, we study the effective susceptibility of second- 
harmonic generation (SHG) in colloidal crystals, which are made of graded metallodielectric nanoparticles 
with an intrinsic SHG susceptibility suspended in a host liquid. We find a large enhancement and redshift 
of SHG responses, which arises from the periodic structure, local field effects and gradation in the metallic 
cores. The optimization of the Ewald-Kornfeld formulation is also investigated. 
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I. INTRODUCTION 



Nonlinear composite materials with high nonlinear optical susceptibilities or optimal figure of 
merit (FOM) have drawn considerable attention for their potential applications, e.g., in bistable 
switches, optical correlators, and so on 1-6 . With several advancements in nanotechnology, such as 
templated sedimentation and dielectrophoresis 7-10 , it is possible to fabricate particles with specific 
geometry. It has also been reported 4 ' 61112 that graded composite materials whose physical prop- 
erties vary gradually in space 13-16 can exhibit enhanced nonlinear optical responses 6 and optimal 
dielectric ac responses 17 , as well as conductive responses 18 . 

It is known that materials lacking inversion symmetry can exhibit a so-called second order 
nonlinearity 1 . This can give rise to the phenomenon of second harmonic generation (SHG), i.e., 
an input (pump) wave can generate another wave with twice the optical frequency (namely, half 
the wavelength) in the medium. In most cases, the pump wave is delivered in the form of a 
laser beam, and the second-harmonic wave is generated in the form of a beam propagating in a 
similar direction. The physical mechanism behind SHG can be understood as follows. Due to the 
second order nonlinearity, the fundamental (pump) wave generates a nonlinear polarization which 
oscillates with twice the fundamental frequency. According to Maxwell's equations, this nonlinear 
polarization radiates an electromagnetic field with this doubled frequency. Due to phase matching 
issues, the generated second-harmonic field propagates dominantly in the direction of the nonlinear 
polarization wave. The latter also interacts with the fundamental wave, so that the pump wave can 
be attenuated (pump depletion) when the second-harmonic intensity develops. In the mean time, 
the energy is transferred from the pump wave to the second-harmonic wave. The SHG effect, 
like the third-order Kerr-type coefficient, involves the nonlinear susceptibilities of the constituents 
and local field enhancement which arises from the structure of composite materials 19-23 . For 
example, Hui and Stroud studied a dilute suspension of coated particles with the shell having 
a nonlinear susceptibility for SHG 19 , and Fan and Huang designed a class of ferrofluid-based 
soft nonlinear optical materials with enhanced SHG with magnetic-field controllabilities 20 . Both 
SHG and phase transformation behaviors can be detected in many nanocrystals by using classical 
methods like X-ray diffraction or Raman spectroscopy 24,25 . Until now, achieving enhanced SHG 
is a challenge 26,27 . 

Theoretical 28 and experimental 29,30 reports also suggested that spherical particles exhibit a 
rather unexpected and nontrivial SHG due to the broken inversion symmetry at particle surfaces 
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(namely, atoms in the surface occupy positions that lack inversion symmetry), despite their cen- 
tral symmetry which seemingly prohibits second-order nonlinear effects. In colloidal suspensions, 
the SHG response for centro symmetric particles was experimentally reported 29 ' 30 . Most recently, 
SHG was also shown to appear for spherical semiconductor nanocrystals 31 . So far, the SHG arising 
from centrosymmetrical structure has received an extensive attention 28-30 ' 32 ' 33 . 

Colloidal crystals have been widely studied in nanomaterials and have potential applications 
in nanophotonics, chemistry, and biomedicine 34 . For colloidal crystals, the individual colloidal 
nanoparticles should be touching since the lattice parameters of the crystals, a, b and c, should 
satisfy the geometric constraint a 2 + b 2 + c 2 = 16Rq, where Rq denotes the radius of an indi- 
vidual colloidal nanoparticle. However, it is possible to achieve a colloidal crystal without the 
particles' touching if the colloidal nanoparticles are charged and stabilized by electrostatic forces. 
In this work, we shall investigate colloidal crystals with the particles' touching. Owing to recent 
advancements in the fabrication of nanoshells 35 ' 36 , we are allowed to use a dielectric surface layer 
with thickness dona graded metallic core with radius a , in order to activate repulsive or attractive 
forces between the nanoparticles. The dielectric constant of the metallic core should be a radial 
function because of a radial gradation, and that of the surface layer can be the same as that of 
the host liquid, as to be used in this work. The latter is also a crucial requirement because oth- 
erwise multipolar interaction between the metallic cores can become important 37 . In this regard, 
the surface layer contributes to the geometric constraint, rather than the effective optical responses 
of the colloidal crystals. In this work, based on the Ewald-Kornfeld formulation [Eq. ©J 38 and 
the nonlinear differential effective dipole approximation (NDEDA) method [Eq. (@])] 39 ' 40 , we shall 
focus on the possibility of achieving such colloidal crystals with desired SHG signals. 

This paper is organized as follows. In section HH we apply the Ewald-Kornfeld formulation to 
derive the local electric field in three typical structures of colloid crystals, and then perform the 
NDEDA method to extract the effective linear dielectric constant and nonlinear susceptibility for 
SHG. In Section Hill we discuss the optimization of the Ewald-Kornfeld summation. Then, we 
numerically investigate the SHG under different conditions in Section [IV] which is followed by a 
discussion and conclusion in Section |V] 
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II. FORMALISM 



Let us start by considering a graded metallic core with radius ao (Fig. [[])• When we take into 
account quadratic nonlinearities only, the local constitutive relation between the displacement field 
D(r) and electric field E(r) is given by 41-43 

AO") = ^2 e v E i( r ) + Xtjk(r)Ej(r)E k (r), i = x,y, z, (1) 

j jk 

where Di(r) and Ei(r) are the ith component of D(r) and E(r), respectively, and Xijk is the 
nonlinear susceptibility for SHG. Here = e(r)Sij denotes the linear dielectric constant, which 
is assumed for simplicity to be isotropic. Both e(r) and Xijk( r ) are functions of r, as a result 
of the gradation profile along the radius r (Fig. [I]). If a monochromatic external field is applied, 
the nonlinearity in the system will generally generate local potentials and fields at all harmonic 
frequencies. For a finite frequency external electric field along z— axis of the form 

E = E (u)e- lujt + c.c, (2) 

the effective SHG susceptibility X2u> can be extracted by considering the volume average of the 
displacement field at the frequency 2uj in the inhomogeneous medium 19,41-43 . In Eq. d2]), c.c. is 
referred to complex conjugate. The graded metallic core can be built up by adding shells gradually, 
making the dielectric constant e(r)(r < a) a radial function, which is schematically shown in 
Fig. Q] We assume that the dielectric constant of the surface and the linear host liquid is a constant 
for convenience, as mentioned in Section HI At radius r, the inhomogeneous spherical particle 
with e(r) and X2u>{f) can have the same dipole moment effect as the homogenous sphere with e(r) 
and X2u{r). The equivalent dielectric constant e(r) can be expressed as the following differential 
equation obtained from the differential effective dipole approximation method 13,39 , 

de(r) (e(r) - e(r))(e(r) + 2e(r)) 



dr re(r) 
On the other hand, the equivalent susceptibility for SHG X2uj(r) can be written as 40 
dx2u(r) _ (2dt ul (r)/dr de 2 ^(r)/dr 2y(u, r) + y(2u, r) - 3 

X2w\j) - — , _ , \ + - — — — + 



(3) 



dr \2e 2 + e UJ (r) 2e 2 + e 2w (r) r 

x 2 (u, r) + ^z 2 (u, r)]x(2u, r) + ^[7x(u>, r) + 2z(u>, r)]z(u, r)z(2u, r) 



rf{2uj,r)f 2 (uj,r) 



(4) 
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where x(u, r) = £2 ^H1 + ^^ , y („ r ) = 2 (g " (r) 7w- ( 1 ) l?o (r) 7^ (r)) , z(w, - e s(r)(Mr)- e „(r)) 



(r)(e w (r)+2e 2 (r))' D ~ 2 e w (r)(e w (r)+2e 2 (r)) ' r > ~ e u) (r)(e u (r)+2e 2 (r)) ' 

and f(u>,r) = Here the indices c<j and 2u; correspond to basic and second harmonics 

for the nonlinear susceptibility, respectively, and X2w{r) denotes the equivalent SHG susceptibility 
of the whole graded spherical particle with radius r. For convenience, we shall denote Xiwi? = a) 
as X2u m the following. 

The above two differential equations [Eqs. ©-©] can be solved numerically as long as the 
gradation profiles are given. For obtaining the effective dielectric constant of the colloidal crystal, 
we refer to the Maxwell-Garnett approximation 44 

e e - e 2 e-e 2 

P - , „ , (5) 



ae e + (3 - a)e 2 e + 2e 2 ' 
where p denotes the volume fraction of the metallic component [see Eq. ©], and the local field 
factor a represents a± (transverse field cases) and an (longitudinal field cases), respectively. Here 
the longitudinal (or transverse) field case corresponds to the fact that the E field of the incident 
light is parallel (or perpendicular) to the uniaxial anisotropic axis. For a± and a\\, there is a sum 
rule 2a ± + an — 3 45 ' 46 . Next, we shall apply the Ewald-Kornfeld model to compute the local 
field factor a for a tetragonal unit cell that can be viewed as a tetragonal lattice plus a basis of two 
nanoparticles. One of the two nanoparticles is located at the corner of the cell, and the other is at 
the body center. Without loss of generality, we consider three representative lattices (Fig. [2]): the 
bet (body-centered tetragonal), bec (body-centered cubic) and fee (face-centered cubic) lattice. If 
the uniaxial anisotropic axis is directed along the z axis, the lattice constants can be denoted by 
a(= b) = tq~^ along x (y) axis and c = iq along z axis, and the volume of the unit cell V c = £ 3 . 
The lattice parameters satisfy the geometric constraint that a 2 + b 2 + c 2 = 16(a + d) 2 , when we 
take into account the dielectric surface layer with thickness d on the graded metallic core. It is 
easy to obtain the value of q from their intrinsic structures: q = 0.87358, 1.0 and 23 represent the 
bet, bec and fee lattice, respectively. Table U shows the calculated values of a± and a\\ versus q, 
according to Eq. © below. The degree of anisotropy of the periodic lattices is measured by how 
q deviates from unity. Here we assume that the colloidal particles are packed closely together. 
Meanwhile, we obtain a relation between q and the volume fraction p of the metallic component, 

with thickness parameter t = (a + d) / a , t > 1. The lattice vector of the tetragonal lattice is 
given by R = £(q~^lx + q~^my + qnz), where I, m, andn are integers. When an external electric 



field E is applied along the x axis, the induced dipole moment p are perpendicular to the uniaxial 
anisotropic axis. Considering the field contribution from all the other particles in the lattice, the 
local field E L at the lattice point r = can be given as 

^ = pEEw) + &cm\ - f £ e ( G >§! «p (^) + t£. ct) 

3=1 R^O c G^O \ / / V 

where ^ = l-(j-l)/2,Rj = |R-[(j-l)/2](ax+6y+cz)|, and 6(G) = l+exp[i(u+v +w)/w]. 
In Eq. (|7]), £> and C are two coefficients given in 46 , B(r) = erfc(^r) /r 3 + 2?7 / (\Znr 2 ) exp(— rfr 2 )) 
and C(r) = 3erfc (rjr) /r 5 + [Arf / (y/7cr 2 ) +6r/ / (y/Tir' 1 )} exp(— rfr 2 ), where erfc(^r) is the comple- 
mentary error function and r] is an adjustable parameter making the summation converge rapidly. 
For details, please see SectionHnl In Eq. ©, p denotes the strength of the induced dipole moment, 
and G the reciprocal lattice vector of R. Thus the local field factor in transverse fields can be 
defined as 

a± = • (8) 

47T p 

For the bet, bec and fee lattices, we obtain a± = 0.95351, 1.0 and 1.0, respectively. Following 
refs 19 and 20, the effective SHG susceptibility of the whole system X2w is given by 

X2» = pX2.T{2u)T 2 {uj), (9) 

where T(u) denotes the factor in a linear system which, for consistency with Eq. © in getting e e , 
should also be determined by using the Maxwell-Garnett approach. Thus, we obtain 

rM = t r , , — • (10) 

V ; (1 - pa)e(u) + (2 + pa)e 2 

Meanwhile, can be obtained through the NDEDA method [Eq. 



III. OPTIMIZATION OF THE EWALD-KORNFELD SUMMATION [EQ. ©] 

In Eq. ©, we see the adjustable parameter r] dominates the accuracy and the efficiency of 
the Ewald-Kornfeld summation 46 ' 47 . An important aspect of the Ewald-Kornfeld summation is the 
tuning in the sense of speed at well controlled errors. It should be chosen carefully in order to make 
the summations both in the real space and the reciprocal lattices converge rapidly 46-48 . On the 
other hand, the r-space cutoff r c and the /c-space cutoff k c are also difficult to be determined 46-48 . 
In this section, we shall analyze the role of these 3 parameters, especially on lattices depicted in 
Fig.[2l Similarly, it can also be applied to other lattice models. 
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In condensed matter physics, complex dielectrics should obey the famous sum rule (3 X + (3 y + 
(3 Z = An, where (3 X , @ y , and f3 z are the local field factors along x, y, and z axes, respectively 
[Eq. (fTTI) l. We shall use this rule to estimate the accuracy of our algorithm. 

First of all, let us investigate how the lattice parameters affect the accuracy of calculation. We 
test 4 lattices with {a = 1, b = 1, c = 1,5, 10} and {a = 10, b = 10, c = 10} respectively. The 
cutoffs in the r— and k— spaces are fixed to 3 which means that we take summation over 6 periods 
in each direction. By using the relation mentioned above, we have 



where the subscript % stands for Cartesian directions x, y, and z. The local-field factors a [Eq. ©] 
and (3 [Eq. (fTTI) l have exactly the same concept, the only difference is that there is a proportionality 
constant An/ 3 between (3 and a, as shown in [Eq. (fTTj) "|. In detail, a± = (An/3)f3 x or (An/3)(3 y , 
and ay = (An/3)(3 Z . 

Next, let us plot the result of the Ewald-Kornfeld summation vs splitting parameter in Fig. [3] 
The total local field factor Ptotai = Px + P y + Pz is normalized by An. When {a = b = c = 1}, 
there is a wide platform with r] goes from 1 to 3. In this case, the correctness of our algorithm 
can be guaranteed by choosing r] with any value within the platform region. However, as c in- 
creases, the anisotropic degree becomes strong, and the original platform rapidly shortens. In the 
lattice {a = b = 1, c = 5}, the proper region of 7] only has the width 0.25. Further, in the lattice 
{a = b = 1, c = 10}, the platform even disappears. At that time, the correct result could not be 
achieved if we would not adjust r c and k c . Another isotropic lattice {a = b = c = 10} is also in- 
vestigated. Although the shape of the figure is similar with that of {a = b = c = 1}, the platform 
is remarkably shrunk. 

Then, let us see how the shape of summation region affects the accuracy. For highly anisotropic 
lattices, i.e. {a = b = 1, c = 10}, the method of cubic region summation is not applicable because 
it counts in many source sites far from the field point along the uniaxial direction, but ignores 
many sites near the field point in the other two isotropic direction (see Fig. HI). It may be improved 
by using a spherical region summation. Set /, m, and n to be sum indices in r— space, and u, v, 
and w to be in fc-space. For the cubic region summation, we just require all these indices to be 
within [-maximum, maximum]. But for the spherical region summation, the requirement turns out 
to be 
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2 Pi ~ 3 
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{alf + {bmf + {cnf < R 2 r , 



(12) 
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{u/af + {v/bf + {w/cf < Rl (13) 

Figure [5] shows the effect of such spherical region summation. The circle line is obtained by 
imposing the restriction that all summation indices are within [-5,5], and the star line is obtained 
by using the spherical region summation. The maximum radii, R r [Eq. (fT2l) l and Rk [Eq. (fT3l) l. 
are the length of the diagonal of the unit cell in r— and A;— space, respectively, 

(al) 2 + (bm) 2 + (cn) 2 <a 2 + b 2 + c 2 = R 2 , (14) 

{u/a) 2 + (v/b) 2 + (w/c) 2 < (1/a) 2 + (1/b) 2 + (1/c) 2 = R 2 k . (15) 

This summation region is smaller than the cubic region, and less sites are evaluated naturally. 
Nevertheless, a wider platform is achieved for this region. The main reason is that the spherical 
region summation tends to count in the sites which contribute much to the field point. The two 
lines in Fig. [5] show the same value when 77 becomes larger. Later we shall see that this is due to 
the result in the /c-space summation that is not modified too much. 

Last, we demonstrate how r c and k c affect the accuracy. In general, the larger radii the cutoff 
has, the more precise the result is. Unfortunately, larger summation regions cost longer computa- 
tion time. Thus, we should take optimized values of r c and k c . In Fig.[6l we calculate the total local 
field factor for four configurations of cutoffs. The left-top graph is copied from Fig.[5l Meanwhile, 
the right-bottom plot is obtained with twice value of r c and k c . We find that by increasing r c , the 
platform would extend to the left with the limit of zero, but there is no such limit by increasing k c . 
As k c increases, the platform could extend to the wide space on the right. 

In conclusion, in order to optimize the Ewald-Kornfeld summation [Eq. ©], we had better 
perform the sum in a spherical shape. Enlarging the summation radius in A;-space is more efficient 
than that in r-space. Our method might test the sum rule first with a large enough r c and k c (here 
time is not the main concern), and then it is convenient for one to choose the center value of r\ 
in the platform. This guarantees the correctness of all the computations performed for achieving 
Table H 



IV. NUMERICAL RESULTS 



For numerical calculations, we set the linear dielectric constant of the nonlinear metallic core 
to have the following Drude form 

e(r) = l-^f- y (16) 



where uj p {t) means a position-dependent plasma frequency, and 7 relaxation rate. For achieving 
the position-dependent plasma frequency, one possible way is to fabricate metallic spherical parti- 
cles containing multilayers each of which is made of different metals. For numerical calculations, 
we take a model plasma-frequency gradation profile 



where is a parameter adjusting the gradation profile. For focusing on the enhancement of 
the SHG response, we take the intrinsic nonlinear SHG susceptibility Xi t0 be a frequency-and- 
position independent real positive constant. 

Figure [7J shows (a) the imaginary part of the effective linear dielectric constant (namely, optical 
absorption), (b) the real and (c) imaginary parts of the effective SHG susceptibility, (d) the modulus 
of X2u/ Xi> an d (e) the FOM (figure of merit) as a function of frequency for longitudinal field cases. 
As C u increases, to p (r) takes on a broader range of value and leads to a broad plasmon band, while 
the plasmon peak shifts to lower frequencies (namely, redshift). The susceptibility of the SHG 
also shows an enhancement and the peak of enhancement can be shifted to lower frequencies, too. 
Generally, the SHG susceptibility and FOM can be enhanced in some frequency regions as 
increases. 

Figure [8] shows the effective responses and FOM as a function of thickness parameter t, for bet 
lattices. For the given lattice it is evident that the effective linear and nonlinear optical responses 
strongly depend on the thickness parameter. Both the redshift and strength of the plasma resonant 
peak are largest at the smallest t for linear optical absorption. In this case, the plasma resonant 
band is also largest for smallest t. Similar behavior can also be found for the nonlinear SHG 
responses. On the other hand, t has an effect on the FOM, too. All of these results come from the 
combination of gradation, local fields, and periodic lattice effects. In fact, the volume fraction for 
different colloidal lattices also contributes to the nature of plasma resonant. For the fee lattice, its 
redshift may lie between bet and bec lattices, and for the current parameters in use, its deference 
between fee and bec (or bet) is smaller than 4% (no pictures shown here). It can also be found that 
as t increases (or, alternatively the colloid crystals become more dilute), the behavior of plasma 
resonant becomes more similar. This can be explained from the dilute limit approximation model 
for the enhancement of optical susceptibility. 

In Figs. |7] and [8] the quantities that can be both positive and negative are plotted in a logarithm 
of modulus. When the quantities pass through zero, the logarithm is very large, thus yielding 




(17) 
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spikes. In addition, we can reach the conclusion that the FOM in the high frequency region is still 
attractive due to the presence of weak optical absorption. 

V. DISCUSSION AND CONCLUSION 

Our main idea is to first reduce the graded metallic cores to effective ones and then consider 
colloidal crystals consisting of such effective nanoparticles embedded in a host liquid that has 
the same dielectric constant as the dielectric shells of the nanoparticles. In doing so , multipolar 
interaction between the metallic cores can become unimportant for arbitrary field polarizations. It 
should be remarked that, since the nonlinear response will depend on the local fields and nonlinear 
susceptibility tensors in the whole structure (due to the broken symmetry at the surface of the 
metallic core), the response can be tensorial and variant within the core. In this work we have 
treated the quantities as a scalar and constant, in order to focus on the effects of lattices and 
gradation of our interest. 

As the value of q increases, the responses in transverse field will have slight difference from the 
longitudinal case (no figures shown here), because of the little difference between a± and a\\ for 
a given structure. From bet, bec to fee lattices, with the increase of q , the longitudinal local field 
factor ay decreases from 1.09299 at bet lattices to 1.0 at bec and fee lattices, while the volume 
fraction p decreases to bec, then increases to fee. This trend can also be seen in other bulk samples, 
such as rhombohedral, orthorhombic and hexagonal. 

We have investigated the cases of graded plasma frequencies, by assuming the relaxation rate 7 
to be a constant. In fact, 7 can also be inhomogeneous. For instance, a position-dependent profile 
for the relaxation rate can be achieved experimentally. One possible way may be to fabricate dirty 
metallic spherical particles in which the degree of disorder varies in the radial direction and hence 
leads to a relaxation-rate gradation profile. In case of graded relaxation rates, the nonlinear optical 
responses can also be adjusted by choosing appropriate gradation profiles for relaxation rates 11 . 

Throughout the paper, the host medium is assumed to be isotropic. It is interesting to see what 
will happen if the host is anisotropic, e.g., for a graded-index host 49 . In this case, the gradation 
is also expected to yield desired enhanced SHG. On the other hand, optical switching in graded 
plasmonic crystals via nonlinear pumping was recently reported 50 , which might also be realized 
in graded colloidal crystals proposed in this work. 

In summary, based on the Edward-Kornfeld formulation, we have theoretically exploited a class 
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of nonlinear materials possessing a nonvalishing SHG susceptibility, which are based on colloidal 
crystals of graded metallodielectric nanoparticles. They have been shown to have an enhancement 
and redshift of SHG signals due to the combination of various effects. 
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TABLE I: Values of a± and ay computed at different q. q = 0.87358, 1, and 2 1 / 3 correspond to bet, bee, 
and fee lattices, respectively. 
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a± 




0.87358 0.953506 


1 OQ9QQ 


0.9 


0.971231 
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1 
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0.999345 


1 001^1 
1 . UU 1 J 1 
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l.UU /4j 


ol/3 
2, i 
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i 
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1.00601 


A QQ7QQQ 

VJ.yo /Voo 


1.4 


1.03492 


0.930155 


1.5 


1.08376 


0.832478 


1.6 


1.15032 


0.699352 


1.7 


1.23137 


0.537268 


1.8 


1.32368 


0.352638 


1.9 


1.42459 


0.150817 
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Figure captions 

Fig. 1. Schematic graph showing the graded metallic core with radius ao embedded in a linear 
host liquid. The metallic core has a dielectric shell that has the same dielectric constant as the 
host liquid, and the core can be built up by adding shells gradually. e 2 denotes the linear dielectric 
constant of the host liquid and shell, and e(r) is the radius-dependent dielectric constant of the 
graded metallic core. 

Fig. 2. Schematic graph showing unit cells of (a) bet, (b) bec, and (c) fee lattices with lattice 
constants a, b, and c, which satisfy a(= b) = £g~ 1//2 and c = Eq. Here, q = 0.87358 (bet), 1.0 
(bec), and 2 1 / 3 (fee). 

Fig. 3. Normalized total local field factors versus splitting parameters in different lattices. The 
rj corresponding to Aotai/(47r) = 1 leads to accurate result of the Ewald-Kornfeld summation 
[Eq. ©]. 

Fig. 4. Cubic and spherical summation region for lattices with different anisotropic degree. 

Fig. 5. Normalized total local field factors versus splitting parameters in cubic and spherical 
region summation. The platform for the i] corresponding to /3 t otai/(4/7r) = 1 is widened by using 
spherical region summation. 

Fig. 6. Normalized total local field factors versus splitting parameters for various r c and k c cut- 
offs. Increasing r c or k c yields the left or right extension of the platform of the r] that corresponds 
to Aotai/(4vr) = 1, respectively. 

Fig. 7. For the bet lattice (longitudinal field), (a) the linear optical absorption Im[e e (u;)], (b) 
Im[X2a,/XiL (c) Re[x2w/Xi]» ( d ) modulus of Xiu/Xi* and ( e ) the F OM=|x 2 J/[XiIm(e e )] versus 
the normalized incident angular frequency of lu/cu p (0) for the dielectric function gradation profile 
[Eq. (fT6l) 1 with various plasma-frequency gradation profiles [Eq. (fT71) 1: C w = 0.3, 0.5, and 0.7. 
Here | • • • | denotes the absolute value or modulus of • • • . Parameters: 7 = 0.02u;p(0), t = 3, and 
e 2 = 2.25. 

Fig. 8. Same as Fig. [71 but for different thickness t = 1.2, 2.0, and 3.0. Parameters: 7 = 

0.02w p (0), C w = 0.5, and e 2 = 2.25. 
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